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We discuss models of computing that are beyond classical. The primary motivation is to unearth 
the cause of nonclassical advantages in computation. Completeness results from computational 
complexity theory lead to the identification of very disparate problems, and offer a kaleidoscopic 
view into the realm of quantum enhancements in computation. Emphasis is placed on the 'power 
of one qubit' model, and the boundary between quantum and classical correlations as delineated by 
quantum discord. A recent result by Eastin on the role of this boundary in the efficient classical 
simulation of quantum computation is discussed. Perceived drawbacks in the interpretation of 
quantum discord as a relevant certificate of quantum enhancements are addressed. 
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I. INTRODUCTION 

There is no known algorithm that may be implemented 
on a classical computer that can estimate the trace of a 
random 2" x 2" unitary matrix efficiently. The scaling 
with n of the number of steps taken by the best known 
classical algorithm to find the trace is exponential. How- 
ever, there is an algorithm that works on a very restricted 
version of a quantum computer [1] that can estimate the 
normalized trace up to a fixed accuracy in a number of 
steps that is independent of n. The specialized quantum 
processor that can run the algorithm has precious little 
that may be identified as "quantum" in the traditional 
sense. Out of the n + 1 qubits used, only one needs any 
amount of purity in its state. All the remaining qubits are 
in the completely mixed state. There is no entanglement 
between the qubit with purity and the rest at any step in 
the algorithm. The "power of one qubit" model, as this 
processor is known, is in the complexity class DQCl [2] 
and is by no means a universal quantum computer. In 
this paper, the DQCl model is of interest because it is one 
of the simplest models in which mixed quantum states 
lead to an exponential improvement over a classical com- 
putation. The aim here is to build a convincing case 
for considering non-classical correlations other than en- 
tanglement as an enabler of computational paths that a 
classical computer cannot take leading to a solution for 
some computational problems in an efficient manner. 

Quantum entanglement is also as a kind of non- 
classical correlation; but not the only kind. It has pro- 
vided deep insights into the role of non-classicality as 
a resource [.>-•"•] in quantum information processing. In 
particular, Josza and Linden [{>] showed that the pres- 
ence of multi-partite entanglement involving a number 
of parties increasing unboundedly with input size is nec- 



essary for a quantum algorithm operating on pure quan- 
tum states to offer an exponential speed-up over classical 
computation. If the entanglement does not increase un- 
boundedly in a computational process, then it can be ef- 
ficiently simulated classically using a classical computer 
to within any prescribed tolerance. An explicit construc- 
tion for such a simulation scheme is also known [7]. On 
the other hand, quantum computations involving a re- 
stricted set of highly entangled states, such as stabilizer 
states, can be simulated efficiently classically, as shown 
by the Gottesman-Knill theorem [8] . Thus entanglement 
appears to be necessary but not sufficient for exponen- 
tial speedups in quantum algorithms implemented using 
pure quantum states. 
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The fundamental problem that we endeavour to ad- 
dress in this paper is of identifying the resource or re- 
sources whose presence, individually or in groups, in 
the quantum states of a quantum information proces- 
sor that enable it to perform tasks exponentially faster 
than the best known classical algorithm running on a 
classical information processing device. We will not be 
looking at sub-exponential speedups provided by quan- 
tum algorithms, as is the case with Grover's search al- 
gorithm [9] and the Bernstein- Vazirani algorithm [111], 
highly interesting as they are. The fundamental ques- 
tion raised above is by no means answered in this paper. 
Instead, the paper is meant to be a survey of what is 
known, and aimed to be an attempt to throw some light 
on how the problem may be approached. Of particular 
interest is the dismantling of the rather artificial divi- 
sion between quantum information processors that work 
using pure quantum states and those that work using 
mixed quantum states, and to identify potential candi- 
dates, unlike entanglement, that may be treated as uni- 
form resources behind the exponential speed ups (if any) 
provided by both types. 
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II. NOT QUITE UNIVERSAL QUANTUM 
COMPUTERS 



While a result along the lines of Josza and Linden prov- 
ing that a quantum resource is necessary and may be even 
sufficient for mixed state quantum computation is very 
much desirable, such a result has been quite elusive as 
well. To simplify the problem, we consider information 
processors that arc not universal, yet perform particular 
computational tasks exponentially faster than the best 
known classical analogue. These quantum computers do 
not usually satisfy all the requirements listed by DiVin- 
cenzo [If] as necessary for assembling a scalable, univer- 
sal quantum computer. They, however, might be consid- 
erably easier to implement than the full-scale, pure-state 
quantum computer discussed by Jozsa and Linden and 
hence they might see practical realization sooner. In the 
remainder of this section, we discuss a few examples of 
such limited quantum information processors which can 
work with mixed or pure quantum states to find solutions 
to specific problems efficiently. 

Suppose we are allowed n identical, non-interacting 
bosons to use in our computation. In particular, we have 
access to only single-photon input states, and photon- 
number-resolving measurements. This is certainly a sub- 
set of the requirements of DiVinccnzo, but this model, 
called boson-sampling, is able to estimate the permanent 
of a n X n unitary matrix efficiently [12]. The perma- 
nent, like the determinant, is a polynomial in the en- 
tries of the matrix. Evaluating the permanent exactly 
is =i^P-complete, a complexity class which can crudely 
be thought of as counting the number of solutions of 
a NP-complete problem [13]. It is possible to estimate 
the permanent of a matrix with non-negative entries effi- 
ciently using Monte-Carlo Markov chains efficiently clas- 
sically [14], but a unitary matrix is beyond the scope of 
these algorithms. 

Another system of limited capabilities is a collection 
of quantum spins whose only mode of interacting is via 
their angular momentum [l-'5]. It has been shown to be 
able to approximate the irreducible representations of the 
symmetric group. There is no known classical algorithm 
for this problem, which is also related to the partition 
function of the Ponzano-Reggc model, a 3-dimensional 
topological quantum field theory on a class of triangu- 
lated manifolds. This computational model, called per- 
mutational quantum computation, is particularly good 
for fault tolerance, as all operations are discrete. In fact, 
the set of n -|- 1 spins interacting via their angular mo- 
menta has only C„ coupling schemes, where C„ is the n*'^ 
Catalan number. Thus, the number of possible two-spin 
interactions is {n + 1)1 C^. Since 
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Stirling's approximation gives the total number of uni- 
taries to be ~ 2°("'°s"). It is known that we need a 
number of unitaries doubly exponential in the number of 



particles involved for the universality of quantum com- 
puting. Thus, permutational quantum computation is 
likely to be less powerful than BQP\ and is not known 
to include all of P either. The final limitation of this 
model, a practical one, is the inability to perform inter- 
fcrometric measurements to extract phase information. 
In spite of the limited power of the permutational model, 
it can approximate the partition function of a physical 
model for which there is no known efficient classical al- 
gorithm [15]. 

A third model of a specialized quantum information 
processor that is of particular interest to us is the 'power 
of one qubit' model. Here the normalized trace of a ran- 
dom unitary matrix is evaluated to any desired accuracy 
in a fixed number of steps independent of the size of 
the unitary. The computational process represented as a 
quantum circuit is given below: 
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The motivation for constructing this model comes from 
quantum computing using liquid-state nuclear magnetic 
resonance (NMR), where it is hard to prepare a number 
of pure qubits in a scalable manner, but easy to get a 
large numbers of them in a very mixed state. 

The quantum state at the end of this circuit, just prior 
to the measurement, is given by 
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wherefrom it directly follows that X and Y measurements 
on the top qubit provide an estimate of the real and 
imaginary part of the normalized trace of the unitary, 
Tr([/„)/2". As the measurement is on a single qubit, 
each individual outcome is either +1 or —1. An average 
over several trials provides the estimate. After L tries, 
the estimate has an accuracy ~ 1/ -s/Z, given by the cen- 
tral limit theorem. The crucial point is that the error 
in the estimate is independent of the size of the prob- 
lem, and this is what makes the estimation algorithm 
efficient. The algorithm can also be thought of as an av- 
eraged phase estimation algorithm. If the input to the 



^ BQP: Bounded error quantum polynomial - the class of problems 
that can be efficiently solved using a universal quantum com- 
puter. P is the class of problems that can be efficiently solved 
using a classical computer 
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lower n qubits is the eigenvector of (7„ with eigenvalue 
e^^^ , the outcome of X and Y measurements on the top 
qubit provides an estimate of cos and sin. cj)j respec- 
tively. The completely mixed state is an equal mixture 
of all the eigenvectors, and the linearity of quantum me- 
chanics immediately leads to the normalized trace, the 
sum of the eigenvalues. 

The top qubit can be replaced by a mixed state given 
by, 

a^^{l2 + aZ). (3) 

X and Y measurements lead to an estimate of 
aTr([/„)/2". This imposes an overhead of l/a^ for 
achieving the same precision in the estimate of the nor- 
malized trace as before. The overhead is tractable as long 
as a is not exponentially small in n. Therefore, we can 
get an exponential quantum speedup with the slightest 
bit of purity in the system. The model can thus be rela- 
beled as the power of the tiniest fraction of a qubit [IG]. 

The power of one qubit model falls into the computa- 
tional complexity class DQCl. The final state of the n + 1 
qubits in this model cannot be represented efficiently us- 
ing the matrix product state formalism for mixed states^ 
in an efficient manner, thus excluding the possibility of 
classically efficient simulations of DQCl computations us- 
ing a large class of algorithms [18]. The power of one 
qubit model has been exploited in a variety of problems, 
from estimating the density of states to the decay of fi- 
delity in chaotic systems, with exponentially enhanced 
efficiency than the best known classical algorithms [1!J- 
24]. 



III. THE COMPLEXITY CLASS DQCl 

The complexity class DQCl represents the family of 
problems that may be solved efficiently with a restricted 
quantum information processor, the 'power of one qubit' 
model. This complexity class is known to be less powerful 
than BQP [2-5], and believed to be more powerful than P. 
The relative inclusions of the different complexity classes 
is shown in Fig. (1). There are several problems that 
have been identified to be in DQCl. We review a few of 
them before addressing the problem of the resources that 
make solving problems in this class using a non-universal 
quantum computer efficient. 

The complexity of estimating the Jones polynomial 
of the plat closure of a braid at the fifth root of unity 
was shown to be BQP-complete [2(i]. The same problem 
for the trace closure of a braid was shown to be DQCl- 




FIG. 1. The relative hierarchy of the complexity classes P, 
BQP and DQCl. The NCI class represents problems that can 
be solved efficiently using a classical parallel computer. 

complete [27], as was the estimation of its HOMFLY'^ 
polynomial [28]. It has been known for quite some time 
that evaluating the partition functions of certain spin 
models, like the Ising model, is related to knot invariants, 
in particular, the Jones polynomial [31]. The difference 
between trace and plat closure is in the way the ends of 
the braids are joined to make a knot. This has a coun- 
terpart in the physical world. Trace closures correspond 
to periodic boundary conditions, while plat closures cor- 
respond to non-periodic boundary conditions [.H2]. The 
partition function of a system depends on its boundary 
conditions. Therefore, one would expect that its esti- 
mation would behave analogously. For the Ising model 
defined on a lattice with periodic boundary conditions, 
estimating its partition function with polynomial accu- 
racy was shown to be DQCl-hard [33, 34]. Of course, 
estimation means approximation to its true value with 
a reasonable accuracy. The notion of 'reasonable accu- 
racy' is a fairly technical one [3o]. As the onc-qubit model 
provides an efficient algorithm for an approximate solu- 
tion to the estimation the Jones polynomial of the trace 
closure of a braid, it should be possible to approximate 
partition functions through the one-qubit algorithm and 
uncover physical reasons for the power of DQCl model. 

The complexity of computing matrix functions like 
traces, determinants, and permanents is an intriguing 
one. We have already seen that the estimation of the 
normalized trace of a unitary matrix is in DQCl. The 
permanent is believed to be extremely hard to evaluate, 
while the determinant of a matrix is easy to evaluate. Al- 
most all the tractable counting problems are, in one way 
or another, reducible to the determinant of a polynomi- 
ally sized matrix. In fact, it is always possible to write the 
permanent of one matrix as the determinant of another. 
The problem is that the size of the latter cannot be poly- 
nomially bounded in the size of the former. A similar 



^ This algorithm [ ] essentially vectorizes the density matrix, and 
looks for an efficient representation for that vector, given by a 
limited Schmidt rank. 



^ The HOMFLY polynomial is a 2-variable generalization of the 
Jones polynomial [ '.']. The estimation of the HOMFLY polyno- 
mial of the plat closure of a braid was also shown to be BQP- 
complete [ ]. 
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dichotomy exists between the trace and plat closures of 
braids, the former of which is DQCl-complete, while the 
latter is BQP-complete. As with matrices, it is always 
possible to transform the plat closure of a braid to the 
trace closure of another, but the error in the estimation 
depends exponentially on the number of strands in this 
case [2(i]. This degrades the quality of the approximation 
of the trace closure exponentially. 

Going beyond the physical and algebraic problems 
mentioned above, we now discuss a combinatorial prob- 
lem, somewhat esoteric, which is in DQCl and reveals 
connections to classical coding theory, graph theory, and 
through that, back to knot theory and statistical physics. 
The problem is that of evaluating quadratically signed 
weight enumerators (QSWEs) [.^fi]. A general quadrati- 
cally signed weight enumerator is of the form 

S{A,B,x,y)= ^ (4) 

b:Ab=0 

where A and B are be matrices over Z2 with B of di- 
mension n X n and A of dimension m x n. The variable 
b in the summand ranges over column vectors in Z2 of 
dimension n, 6"^ denotes the transpose of b, \b\ is the 
weight of b (the number of ones in the vector 6), and all 
calculations involving A, B and b are modulo 2. The ab- 
solute value of <S'(^, B, x, y) is bounded by {\x\ + |?/|)". In 
general, one can consider the computational problem of 
evaluating the sum in Eq. (4). In particular, it is known 
that for integers k,l, evaluating S{A, B,k,l) exactly is 
#P complete. Now, let A be square, of size n x n, and 
let A denote the lower triangular part of A, which is the 
matrix obtained from A by setting to zero all the entries 
on or above the diagonal. Let diag{A) denote the diago- 
nal matrix whose diagonal is the same as that of A and / 
denote the identity matrix. For matrices C and D with 
the same number of columns, let [C;D] denote the ma- 
trix obtained by placing C above D. Then, the following 
two theorems hold [3G] 

Theorem 1 Given diag{AJ = /, k, I positive integers, 
and the promise that \S{A, A,k,l)\> (fc2+?2yV2/2, de- 
termining the sign of S{A, A,k,l) is BQP-complete. 

Theorem 2 Given diag[A) = /, k,l positive inte- 
gers, and the promise that \S{[A; A^], A,k,l)\ > {k^ + 
Z2)»/2/2^ the sign of S{[A; A^],^ , k,l) can be efficiently 
determined by the the power of one qubit model. 

In these two problems, the integers k and I can be re- 
stricted to 4 and 3, respectively, without affecting their 
hardness with respect to polynomial reductions (using 
classical deterministic algorithms) . Note that in the sec- 
ond one, the question of DQCl completeness is still open. 

Lidar made the first connection between QSWEs and 
partition functions of spin models on certain classes of 
graphs [37]. In that case, A denotes the adjacency ma- 
trix of the relevant graph. Though these classes are 



highly restricted at present, attempts at broadening their 
scope are ongoing. Through the connection to QSWEs, 
DQCl can be applied to analyze problems in coding the- 
ory, like weight generating functions of binary codes. 
Studying the complexity of such problems will certainly 
give us a better understanding of the complexity classes 
P, BQP, and DQCl and their relative inclusions. As 
shown in Fig. (2), the 'power of one qubit' model can 
be employed to tackle problems in graph and knot the- 
ory, disciplines that have provided scores of problems 
for theoretical computer science in general and quan- 
tum computation in particular. Identifying problems in 
DQCl and studying them in relation to BQP-, P-and NP- 
complete graph and knot theoretic problems might un- 
ravel the underlying mathematical structures that sepa- 
rate P, BQP, and DQCl. Such an understanding should 
illuminate greatly the mathematical framework underly- 
ing the power of quantum computation, just as study- 
ing the complexity of partition functions might unravel 
the physical foundations for that power. The diagram in 
Fig. (2) outlines prospective avenues for exploring such 
problems in the context of mixed-state quantum compu- 
tation, and has the potential of opening up completely 
new avenues in quantum computation research. 



IV. NON-CLASSICAL CORRELATIONS IN 
DQCl 

In the' power of one qubit' model discussed earlier, the 
top qubit is always separable from the rest [!()]. This lies 
at the centre of all the discourse involving the 'power of 
the one qubit' model. If entanglement is the sole form 
of quantum correlations capable of providing quantum 
enhancements in computation, and there is none of it 
between the top qubit and the rest, how does the infor- 
mation travel from the bottom qubits to the top, allowing 
us estimate to the trace of the unitary that acts on the 
bottom qubits alone? Looking for other quantifiers of 
non-classical correlations in a quantum state is therefor 
a good starting point for addressing this problem. 

Quantum discord was proposed as a measure of purely 
quantum correlations [38, 39] about a decade ago. It was 
presented mathematically as a difference of two classi- 
cally equivalent, but quantum mechanically inequivalent 
forms of the mutual information. The first is the quan- 
tum mutual information, which is a measure of all cor- 
relations - quantum as well as classical - in a quantum 
state. For two systems A and B, the mutual information 
is 

I{A: B) = H{A) + H{B) - H{A,B). (5) 

If A and B are classical systems whose state is described 
by a probability distribution p{A,B), then H{-) denotes 
the Shannon entropy, H{p) = —J^jPj^'^SPj^ where p 
is a probability vector. If A and B are quantum sys- 
tems described by a combined density matrix pab, then 
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FIG. 2. Problems that have connections to the DQCl complexity class: Black lines show the connections already known in the 
literature. The blue lines show some of the possible connections that can be made. 'Trace' and 'Plat' refer to closures of braids 
whose Jones polynomials are of interest, while 'Periodic' and 'Aperiodic' refer to boundary conditions for Ising models. 'Graph 
enumeration' refers to the class of graph theoretic problems like the evaluation of the Tutte, the dichromatic, or the chromatic 
polynomial. The Kauffman bracket polynomial is a knot invariant which is a generalization of the Jones polynomial, and is 
related to partition functions of more esoteric statistical mechanical models like the vertex model and the IRF (interactions 
round a face) model. 



H{-) stands for the corresponding von Neumann entropy, 
Hip) = -Tv{p\ogp). 

For classical probability distributions, Bayes's rule 
leads to an equivalent expression for the mutual informa- 
tion, I{A : B) ^ H{B) - H{B\A), where the conditional 
entropy H{B\A) is an average of Shannon entropies for 
i3, conditioned on the alternatives for A. For quantum 
systems, we can regard this form for I{A : B) as defin- 
ing a conditional entropy, but it is not an average of 
von Neumann entropies and is not necessarily nonneg- 
ative. An alternative route to generalizing the classical 
conditional entropy to the quantum case is to recognize 



that classically H{B\A) quantifies the ignorance about 
the system B that remains if we make measurements to 
determine A. When A is a quantum system, the amount 
of information we can extract about it depends on the 
choice of measurement. For a measurement given by 
a positive-operator- valued measurement (POVM) {Ilj} 
corresponding to outcomes j , the state of B after a mea- 
surement on A is given by 

PB\j = TvA{^jPAB^j)/Pj, Pj = Tt a.b{pab^])- (6) 

A quantum analogue of the conditional entropy can then 
be defined as H^Yi.j{B\A) = J2j Pj^iPB\j) > 0- Since 
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Pb = '^jPjPB\ji the concavity of von Neumann entropy 

implies tliat H{B) > H[Yi^y{B\A). We can now define 
another quantum version of the mutual information, 



J{n,}(A : B) EE H{B) - i/{n,}(B|^) > 



0. (7) 



Performing projective measurements onto a complete set 
of orthogonal states of A effectively removes all non- 
classical correlations between B and A. In the post- 
measurement state, mutually orthogonal states of A are 
correlated with at most as many states of B. It is easy 
to see that these sorts of correlations can be present in 
an equivalent classical system. If on the other hand, the 
system B is bigger than A, the correlations are certainly 
non-classical, or quantum. 

The value of J7{nj}(^ ■ B) in Eq. (7) depends on 
the choice of {Hj}. We want Jsji.y{A : B) to quan- 
tify all the classical correlations in pab, so we maximize 
J^[Yij}{A ■ B) over all {11^} and define a measurement- 
independent mutual information J{A : B) = H{B) — 

H{B\A) > 0, where H{B\A) ee min^n,} Ej Pj^(Psb) 
is a measurement-independent conditional information. 
Henderson and Vedral [4U] investigated how J^{A : B) 
quantifies classical correlations, and an operational in- 
terpretation was also provided [41]. The criteria they 
postulated for a measure of classical correlations C were 

1. C = for product states 

2. C is invariant under local unitary transformations. 
This is because any change of basis should not affect 
the correlation between two subsystems. 

3. C is non-increasing under local operations. If the 
two subsystems evolve independently then the cor- 
relation between them cannot increase. 

4. C = H{pa) = H{pb) for pure states. 

One can then define a measure of purely quantum cor- 
relations as the difference of the total correlations in a 
system and C. Taking the quantum mutual information 
to be a measure of total correlations in a system, the 
purely quantum correlations can be measured by 



V{A, B) EE I{A : B) - J {A : B) 
= H{B\A) - H{B\A). 



(8) 



here is quite a substantial literature on the properties 
of quantum discord, its evaluation in a variety of sys- 
tems, and its role on various phenomenon other than 
computation, and has been reviewed in this volume [42]. 
Quantum discord for pure states reduces to quantum en- 
tanglement, given by the von-Neumann entropy of the 
reduced state. Secondly, quantum discord is asymmet- 
ric between the two parties, 'D{A,B) ^ T>{B,A). Fi- 
nally, the maximum value of quantum discord is the von- 
Neumann entropy of the measured subsystem, in this 
case 'D{A,B) < H{A). The set of states with zero dis- 
cord will form an important ingredient in our discussion. 



There are several equivalent characterizations of this set, 
but an helpful one is that a state has zero discord if and 
only if it is block-diagonal in the marginal eigenbasis of 
the measured subsystem [43-46]. 

It was proposed that quantum discord might be the 
resource behind the speedup offered by the power of one 
qubit model [47]. The reason is two- fold. Firstly, there 
is no quantum entanglement between the top qubit and 
the rest, while there is a non-zero amount of quantum 
discord. Secondly, the amount of quantum discord across 
this split is a constant fraction of the maximum possible 
for this system, which is 1 ebit. To better understand the 
role of discord in the DQCl model, let us generalize the 
circuit slightly. The state of the top qubit is as in Eq. (3) 
and it is not pure. The Hadamard gate H is replaced by 
a general 1-qubit unitary V given by 



V 



e"^ cos e 
— e*^ sin( 



e'^ sin 9 
e-^f cose 



(9) 



Then the input state is 

1 / l + acos26' -ae*'' sin 261 
2"+i V -ae"''^sin2e' l-acos26' 



where 6 = ip + a-nd the output state, in block form, is 

Pn+i[a) - -(ae-**sin26')C/t (1 - acos26l)/„ 

(10) 

The circuit looks like 
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and the state is completely separable across the split 
given by the dotted lines. If the unitary has the eigende- 
composition J7„ = Ej ^"^^ kjXsjIj the state can be writ- 
ten in the form 



p„+i(a) = -l^^(ie,)(oi + ic,)(o 



(11) 



where 

l^j) = ajO) + e^(^^+*)611), 10) = cjO) + S'^^+^^d\\), 
with the real parameters a, &, c, d given by 



7 



1 + a cos 26*, 



1 — a cos 29, ab + cd = a sin 26*, 



The state is thus separable across the dotted boxes above. 

However, there is quantum discord across that split for 
a typical unitary [47], when the measurement is made 
on the subsystem A. One can show that the optimal 
measurement lies in the X — Y plane, given by DJIa = 
cos(f>X + sin (j)Y. In Fig. (3), the difference I — J' is 
plotted as a function of the measurement angle </>, for the 
unitary corresponding to the trace closure of a 4-strand 
braid [5(J]. The minimum of the graph, which gives us 
the value of discord is also shown. 




FIG. 3. The quantum discord in the circuit that estimates 
the Jones polynomial of the trace closure 4-strand, with the 
first and the second strands crossing. The discord is D — 
0.23, and we plot the difference of the mutual and conditional 
information as a function of the measurement DJIa ~ cos (j)X + 
sin^y. The optimal measurement is given by = 0.31. 



The modified version of the power of one qubit model 
that we have considered here provides approximations to 
the value the Jones polynomial of a simple knot at the 
fifth root of unity. A recent experiment^ has confirmed 
that there is indeed non-zero discord in the DQCl cir- 
cuit that estimates the Jones polynomial of the trace clo- 
sure of the 4-strand braid, with the first and the second 
strands crossing [50, 53]. In this experiment, the polar- 
ization of the top qubit was quite small, a ^ 10~^, [52] 
which leads to a very small, but strictly non-zero value 
of discord. The plot in Fig. (3) is for a = 1. The pres- 
ence of discord in this circuit together with the absence 



I 

of entanglement has led to the proposal of quantum dis- 
cord being the resource for the enhanced power of the 
class DQCl. However, this is not a formal proof that 
discord is necessary for quantum enhancements in the 
DQCl model. 

As we mentioned earlier in this section, discord is 
asymmetric across the subsystems, and the quantum dis- 
cord in pn+i{a) for measurements made on subsystem 
B, given by the projectors onto the eigenvectors of J7„, 
is zero. This can be used to repudiate the above claim 
of discord being a resource for the power of DQCl. To 
explore this further, let us perform a simple calculation. 
Say we measure the n qubits in _B, using an operator 971b. 
The expectation value of such a measurement is 



Tr 



Tr (Ml 



Mb 
Mb 



(12) 



An all-optical experimental implementation of the DQCl model 
has also been carried out [-"ji], as has been liquid-state NMR 
implementations to evaluate Jones polynomials [ix, I')]. 



The result is independent of the input state a and the 
unitary V. This shows that no information about the sys- 
tem A is obtainable by measuring system B, no matter 
what measurement we make. Since any computation can 
be thought of as the flow of information, zero discord in 
one direction in the one qubit model merely reflects the 
just derived fact that there is no information propagat- 
ing across the dotted boxes from A to B. Thus, what 
might be considered as a drawback for quantum discord, 
its asymmetry, can heuristically be understood as des- 
ignating the direction of propagation of quantum infor- 
mation. Indeed, pure-state quantum computations are 
completely unitary, thereby making information transfer 
between parts of such a system necessarily symmetric. 
For such systems, quantum discord is identical to quan- 
tum entanglement which, as a measure of correlations, 
is also symmetric between the systems involved. Not so 
for a mixed-state quantum computer. As in the power of 
one qubit model, states with non-zero quantum discord 
might provide a possible explanation to quantum compu- 
tation with mixed quantum states, and the asymmetry 
of the measure could be vital to its veracity. 

Recently, it has also been shown that almost all states 
possess non-zero quantum discord [43]. This is unlike 
entanglement, where the measure of the set of separable 
states is non-zero. In particular, there is a ball around 
the completely mixed state inside which all states are sep- 
arable [54, 55]. This has generated inordinate amounts 
of excitement, in its manifestation as the sudden death 
of entanglement. There is no such phenomenon for quan- 
tum discord, and the set of states with zero quantum dis- 
cord is of measure zero. Consequently, one might wonder 
how so ubiquitous a quantity as quantum discord might 
be a resource. 
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One may build some intuition on the role of nonclassi- 
cal correlations in quantum computing using the follow- 
ing heuristic: Consider a generic quantum information 
processor that uses registers of qubits for input and out- 
put, and an additional register of qubits for the computa- 
tion. Both the input and output are human-readable bit 
strings; the input bit string being loaded onto the qubit 
register typically by a process or quantum state prepa- 
ration and the quantum state of the main register at the 
output being reduced to a bit string by a measurement. 
The initial state of the computing qubits may be chosen 
so that the starting state of the quantum computer has 
no nonclassical correlation, including discord, across any 
bipartite or multipartite division. The state at the out- 
put end may also be chosen likewise by measuring them 
if needed. 

One can now imagine the computational process to 
be a sequence of discrete steps that takes the all the 
qubits from the initial, zero-discord, state to the final, 
zero-discord, state. The discrete steps correspond to the 
action of each gate set in the computational sequence. 
Since a comparison with classical computing is attempted 
wherein the computational steps are discrete transforma- 
tions and not unitary gate actions, let us ignore what 
happens while each gate is acting. In other words, we 
consider only the states before and after the action of 
each gate. Imagine that we now run a classical algo- 
rithm as well as a quantum algorithm for the solving the 
same computational problem on the quantum informa- 
tion processor. The classical algorithm will demand that 
at each step the state of the quantum computer have 
no nonclassical correlation across any bipartite division. 
Note that a deterministic classical algorithm will be even 
more restrictive, demanding product states in the com- 
putational basis at each intermediate step. So we are 
including probabilistic classical algorithms as well. This 
means that one has to go from the initial to the final state 
jumping from one zero discord state to another without 
ever leaving this set of measure zero. A true quantum 
computation however, lets one take alternate paths in- 
volving intermediate states with nonclassical correlations 
which might result in reaching the desired output state 
in significantly - perhaps even exponentially - fewer com- 
putational steps. 

This above heuristic is represented in Fig. (4). A classi- 
cal algorithm will restrict the state of the quantum com- 
puter on the blue lines representing zero discord states. 
A quantum algorithm will not. The states pt and pf are 
the initial and target states of the quantum computer. 
Attempts are ongoing to formalize this simple intuition. 
It is to be noted that this is more of a qualitative than 
a quantitative argument, and what appears to be im- 
portant here is the distribution of zero discord, classical 
states in the space of all quantum states. 

Eastin [50] has recently proved a significant result 
along the lines described above. He considered a quan- 
tum computation for which the input state is a product 
state diagonal in the computational basis, p'^ = ®kPk- 




FIG. 4. A schematic representation of the space of quantum 
states. The gray space S is the set of separable states, the rest 
is the set of entangled states E. The larger ellipse is the set 
of all states. The blue lines are the set of zero discord states, 
and they intersect at the completely mixed state. If one is 
proceed from pi to pf touching down only on the blue lines, 
the computation might take more steps than the case where 
the state can explore the entire space. 



The computation is done by a sequence of unitary gates 
{G*}, and concluded by single-subsystem measurements 
determining the outcome of the computation. Eastin 
showed that if the quantum state of the qubits involved 
at every step of the computation is concordant then the 
computation can be simulated efficiently using a clas- 
sical computer. A concordant state is one that is di- 
agonal in a product basis. This is equivalent [4()] to 
ViA^A) = 0_= 'D{A,A) for any subsystem A and its 
complement A. In other words, the quantum discord of 
the state across any bipartition must be zero. 

Eastin's result is restricted in the sense that it docs 
not allow any nonclassical correlations, even those lim- 
ited to a subset of qubits, during the computation. This 
is unlike the Joza-Linden result which showed that the 
computational steps can be efficiently simulated even if 
entanglement, restricted to a subset of qubits, is present. 
Eastin's result is also a qualitative one in which the pres- 
ence and absence of discord is considered and not the 
amount of discord. In addition to this, Eastin defines 
quantum discord as a minimization over a set of com- 
plete, orthogonal set of one dimensional projectors [■)!)], 
and not general POVMS [57]. The gate operations in the 
computation are restricted to one and two qubit gates as 
well. This restriction, if relaxed, could lead the way to a 
more quantitative version of Eastin's result, in that one 
could study the classical simulation of systems which are 
approximately concordant. It would also allow for study- 
ing the classical simulation of systems with quantum dis- 
cord restricted to blocks of bounded size. This would 
be akin to the results of Vidal [7] , as extended to mixed 
states. 

Finally, for completeness, we consider an argument 
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that was raised against considering discord as a resource 
behind exponential spcedups. Consider the DQCl circuit 
when 

Un = kUl (13) 

This is a necessary and sufRcicnt condition for vanishing 
discord in the pre-measurement state in the DQCl model. 
Let us first note that k = e'^^ , must be a phase. It 
has been suggested that for a unitary of this form, the 
DQCl circuit provides an exponential speedup in the ab- 
sence of quantum discord, and so quantum discord can- 
not be the source for the quantum advantages in this 
model [53]. However, if we decompose the unitary C/„ in 
terms of any universal gate set, the above condition does 
not impose any restriction on the intermediate states of 
the computation. Therefore, quantum discord may very 
well be present in any of the states preceding the final 



one. As a simple example, consider Un = Z, the Pauli 
operator. This operator can be composed of two phase 
gates P, since Z ~ P^, where 

.^(;°). (») 

Although, the final DQCl state has no quantum discord, 
the intermediate state after the application of one phase 
gate has non-zero quantum discord. Thus, zero quantum 
discord in the final state does not imply the absence of 
discord at all stages of the computation. Furthermore, 
zero quantum discord in every intermediate state does 
not imply the same at the final state. This is because 
the unitaries satisfying Eq. (13) do not form a group, 
implying that if we imposed this condition at each step of 
the decomposition, the final unitary need not, in general, 
have this property. In other words, if Un = W1W2 ■ ■ ■ Wt, 
and each Wi satisfies Eq. (13) (upto the phase factor fc). 



Un = WiW2---Wt 



wlwl ■ 



■Wl ^ {WiW2---Wt)^ ^Ul 



(15) 



It can thus be concluded that a DQCl-circuit with a uni- 
tary satisfying Eq. (13) is not a valid rebuttal of the 
suggestion that quantum discord might be responsible 
for the speedup in the power of one qubit model. 



V. CONCLUSIONS 

Identifying the boundary between quantum states with 
nonclassical correlations and those without would signifi- 
cantly aid in designing information processing algorithms 
and devices that can provide exponential speedups over 
their classical analogues. We have reviewed the role of 
quantum discord as a quantifier of nonclassical correla- 
tions. We have presented some evidence that discord 
might be necessary for quantum speedups in mixed-state 
quantum computations. In particular, we noted that the 
formal proofs in this direction are more qualitative in na- 
ture than quantitative, and a direct connection between 
the amount of quantum discord in a system and its com- 
putational power is lacking. Parallel to developing the 
connection between quantum discord and computation, 
there have been advances in understanding the nature of 
quantum discord as well. Quantum discord has an op- 
erational interpretation in quantum information theory 
in terms of quantum state merging [58, 59], as well as in 
terms of the thermodynamics of quantum and classical 
Maxwell's demons [GO, iil]. 

There are results [15] that suggest that quantum dis- 
cord is not sufficient for quantum speedups using mixed 
states. The DQCl version of permutational quantum 
computation provides an estimate of the normalized char- 
acter of the irreducible representations of the symmetric 



group. This problem is known to be in BPP, since a ran- 
domized sampling algorithm can provide the answer, yet 
there is quantum discord in the circuit that implements 
it. Additionally, this is indirect evidence that permuta- 
tional quantum computation is weaker than DQCl [15]. 
There are no results yet on the relative strengths amongst 
classes of non-universal quantum computational models, 
and classes of problems like counting, decision, or sam- 
pling. Such results will be extremely valuable in under- 
standing features of quantum mechanics, and their im- 
plications as quantum resources. 

The unlikelihood of a single quantity certifying all 
forms of quantum advantages is becoming apparent with 
our increased understanding of quantum information sci- 
ence. Yet, we can attempt at formalizing the notion of 
quantumness into as few an independent number of quan- 
tities as possible for a given system or scenario. In course 
of this exercise, we must make every attempt to keep our 
system or scenario as generic as possible. What is pro- 
gressively becoming clear is the presence of a trade-off 
between the conciseness of our understanding, and the 
scope of its validity. In the long run, the validity will be 
governed by the practical feasibility of the model. Given 
that the construction of a full-scale quantum computer is 
still a challenging prospect, intermediate models of com- 
putation are more than theoretical curiosities. If such 
models are able to provide (even polynomial) speedups 
over best known classical algorithms for certain tasks, it 
would be a significant advance in the field. More so if the 
tasks are widely relevant, for instance, matrix inversions, 
root-finding, Fourier transforms. Therefore, intermedi- 
ate models of computation are a vital component in the 
quest for tangible quantum advantages. In the process. 
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they will clarify what differentiates quantum from classi- 
cal when it comes to computational advantages. 

Quantum computing with mixed states seem to be 
a general enough setting for which we ought to find 
the certificates of quantum enhancement. This is more 
than a mere theoretical exercise with abstract implica- 
tions. Substantial experimental effort in quantum com- 
putation is directed at maintaining the purity of the 
quantum systems. A better understanding of mixed- 
state quantum computations that provide exponential 
enhancements will help direct valuable resources in har- 
nessing the available potential, rather than shying away 
from them merely because they involve states that are 
not pure. DQCl is such a model. The promise of the 
DQCl model lies in its strong interconnections to several 
areas of physics and mathematics. The set of complete 
problems for this class is varied, from knot theory to cod- 
ing theory, and statistical mechanics to topological field 
theory [62, 63]. This makes it an ideal workbench in the 
analysis and understanding of the power of different phys- 
ical models and the computations of different invariants 
related to these models, as well as potential resources for 
quantum enhancements such as quantum discord. The 
interplay between the computational efforts needed to es- 
timate elements versus characters of irreducible represen- 



tations, individual matrix elements versus traces of uni- 
tary matrices, knot invariants of plat versus trace closures 
of braids, and partition functions of systems with peri- 
odic versus non-periodic boundary conditions on square 
lattices provide us with a variety of possibilities to unlock 
the true sources of quantum enhancements in computa- 
tion, as indicated in Fig. (2). These are all problems 
primed to be explored. 
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